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Abstract 



Leading-order electromagnetic loop corrections to the one-pion-exchange 
potential are computed within the framework of chiral perturbation theory. 
These corrections generate an effective nucleon-nucleon potential, V^ 7 , which 
supplements the sum of OPEP and the nucleon-nucleon Coulomb potential. 
This potential is charge dependent and its construction is demonstrated to be 
gauge invariant. The potential has been included in the Nijmegen partial- 
wave analysis of np data. A particular renormalization scheme is chosen that 
leads to a negligible change in the tt NN coupling constant and in the np 
scattering length and effective range. 
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Isospin violation is an august topic ^ |2|] in nuclear physics, and considerable 
progress has been made in understanding the consequences of isospin violation in 
nuclei, particularly in few-nucleon systems^, |J. Nevertheless, a comparable under- 
standing of the underlying mechanisms for isospin violation in the nuclear force has 
been slower in developing. It is our purpose here to present for the first time complete 
analytic results for one such mechanism obtained using chiral perturbation theory, 
and to determine its effect on the ^NN coupling constant and the np scattering 
length and effective range. 

It has been known for decades that the various mechanisms responsible for isospin 
violation could be graded in strength, with charge- dependent (CD) nuclear forces 
being generally stronger than charge-symmetry-breaking (CSB) nuclear forces, while 
strongest of all is the long-range Coulomb force between protons. Only recently, 
however, have these empirical observations been linked quantitatively to microscopic 
strong-interaction mechanisms based on symmetries |5], ||, [?]]. 

The strength of isospin violation can be understood |J using dimensional power 
counting, a property associated with chiral effective Lagrangians||. This technique 
expresses the strength of any nuclear interaction in terms of several energy scales: a 
heavy scale Mqcd ~ 1 GeV (the characteristic QCD mass) and the lighter scales f n = 
92.4 MeV (the pion decay constant), = 139.6 MeV (the scale of explicit chiral- 
symmetry breaking), and Q (the effective momentum in a nucleus, which can be taken 
to be ~ m n ). In addition, up-down quark-mass-difference-induced isospin violation^ 
carries an extra factor of e = (m^ — m u )/(md + m u ) ~ 0.3, while electromagnetic 
interactions carry powers of the fine-structure constant a (~ 1/137). Dimensionless 
strong-interaction constants of "natural" [§] size ~1 (such as qa = 1-26, the axial- 
vector coupling constant) are ignored. The general technique is not dependent on 
any particular model, but specific force models can be analyzed and expressed in 
terms of these quantities, as well. Signs are not determined by such arguments. 

Using these techniques one finds that the dominant (isospin-conserving) compo- 
nent of the nuclear force in light nuclei, OPEP, has a strength 0] (corresponding to 
an energy shift) ~ 15 MeV/nucleon pair, while the pp Coulomb interaction (CD and 
CSB) is roughly 1 MeV in strength. The small Breit corrections to the latter (viz., 
magnetic and velocity-dependent forces) have a nominal 20 keV strength, as does the 
effect of the nucleon-mass difference. The difference in pion masses in OPEP gen- 
erates a CD nuclear force of roughly 1/2 MeV in strength. Simultaneous exchange 
of a pion and a photon between nucleons has a nominal 30 keV strength, as does 
the electromagnetic modification (at the quark level) |10[ of the tcNN coupling con- 
stant; both are CD mechanisms. The leading-order CSB force generated by the effect 
of differing quark masses on the kNN coupling constants may be somewhat larger 
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~ 90 keV. These estimates (which could easily vary by a factor of two) correspond 
rather well to those obtained from the well-studied isospin violation in few-nucleon 
systems 0, f|, p7o[| . Thus, one finds that CD is larger than CSB, and both of these 
mechanisms are larger than the tiny isospin- violating forces of Type IV Q, all of 
which was first demonstrated in Ref. || using dimensional power counting. 



Isospin violation in the nuclear force has also been studied as a natural extension 



of the partial- wave analysis (PWA) of both pp and np scattering data|TTJ, together 



with the tiny amount of nn scattering information ||12||. It has long been known that 
the AW scattering lengths are all different, and their differences scale roughly as 
the interaction strengths estimated above. These PWAs have allowed the 7c°pp, n np, 
and ^NN coupling constants to be separately determined ^] , and the Goldberger- 
Treiman (GT) discrepancy |14|] to be obtained, as well. At the present time||T3|, |10 | 



both the CD and CSB ttNN coupling constants are consistent with zero (within < 1% 
overall uncertainty in the ttNN constants), and the GT discrepancy (d — 1), which is 
isospin conserving but chiral-symmetry breaking, is approximately 2%. 



The Nijmegen PWA [|TT|1 incorporates pion-mass differences in OPEP (K-), nucleon- 
mass differences, and the full electromagnetic interaction |1| (including vacuum po- 
larization) when it determines the ttNN coupling constants. Only the effective tt-j- 
exchange interaction has been missing, and we have seen that its strength is com- 
parable to that of several isospin-violating mechanisms normally incorporated. The 
nominal strength of that force is (as we shall see) (a/n)V n , and this is comparable to 
the uncertainty in the ^NN coupling constant, f 2 = (gA m w+ ^/2/^) 2 /47r. Therefore, 
the task that remains is to calculate the effective 7r-7-exchange force, incorporate it 
into the Nijmegen PWA, and redetermine f 2 . 

The calculation is performed within the framework of chiral perturbation theory ||, 
^|, [TB| to leading order in 1/Mqcd- This allows us to make the static approximation 



for nucleons (M^r — > oo), and restricts us to the leading-order terms (in a chiral 
expansion) in the effective Lagrangian. Because we are generating EM corrections to 
one-pion exchange between nucleons (short-range forces were discussed previously^], 
|T0"ll), only single- loop corrections involving photons are required. In static order, the 
necessary elements of the leading-order interaction Lagrangian are: 

L (0) = -^iV[cr.V(t-7r)]iV + ^iVo--A(tX7r)3Ar 

Jn Jit 

-eA N{\ + t 3 )N - eA^n x ^tt) 3 + \e 2 A\-K 2 - tt|) + • • • , (1) 
where we have used Cartesian notation for the pion charge states (iTi) but have oth- 



erwise conformed to the notation and conventions of Refs. [17| and JT0[. Additional 
terms in the Lagrangian that are not required here can be found in the latter refer- 
ence. The nucleon (Pauli) spin operator is cr, its isospin operator is t a (t 3 = ±|), e 
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is the fundamental (proton) charge, and N 1 is the photon field. Although to lead- 
ing (static) order the ir^yN seagull interaction involves only A and the nucleon EM 
interaction involves only Aq, the full pion EM interaction is required. 
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Figure 1: Electromagnetic loops contributing to isospin violation in the OPE nuclear 
force in leading order. Solid lines are nucleons, dashed lines are pions, and wavy lines 
are photons. 



The calculation requires all diagrams of Fig. (1). Figure (la) subsumes 4 separate 
diagrams (a bubble on each nucleon leg), (lb) 2 diagrams, (lc) 4 diagrams, (Id) 2 
diagrams, (le) and (If) each 2 diagrams, (lg) is unique, while (lh) is accompanied by 
an EM tadpole (obtained from the last term in Eq. (1)) that is not shown. Altogether, 
19 diagrams were separately calculated in both Coulomb and Feynman gauges. All 
graphs differ in the two gauges. Coulomb gauge has the advantage of eliminating a 
priori any infrared divergences, but considerable complexity is introduced into the 
calculation. Infrared divergences arise in Feynman gauge but cancel, while a compli- 
cated set of terms that ultimately cancel in Coulomb gauge never arise in Feynman 
gauge. The sum of diagrams is the same in both gauges, a natural consequence of 
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any complete calculation, rather than one focusing on a few selected mechanisms. 



In order to obtain physical and useful results we have to use a renormalization 
prescription. We follow Ref. |1| and keep only nonanalytic terms in divergent loops, 
incorporating all analytic terms into the definitions of the coupling constants and 
masses. The hard parts of vertex loops renormalize the CD irNN coupling constant 
(denoted by Pio 



in Refs. 
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for example. In addition, diagrams (le) and (If) 
contain a component corresponding to the iteration of the static OPEP and Coulomb 
potential (in the general case) and this has been subtracted from the full ampli- 
tude, since it is automatically included when solving the Schrodinger equation. The 
remaining amplitude (SV7) is still gauge invariant and the relation 



S. 



7T7 



-iV 7 



7T7 



(2) 



defines the effective 7r-7-exchange potential, V n ^, between two nucleons arbitrarily la- 
belled 1 and 2. Its isospin structure allows only charged-pion exchange and therefore 
affects only np scattering. The usual OPEP corresponding to charged-pion (tt c ) ex- 
change, V 7Tc , is similarly obtained. Their sum is given by the simple momentum-space 
expressions 



K 7 (q) + K c (q) 

V^{(3) 



t x • t 2 - *M)(«ri • q<r 2 • q) [K 7 (/3) + VM] , (3) 



a 

7i 



2\2 



(1-/3 2 ) 



2/3 4 (l + 

where OPEP is determined from 



ln(l + /3 2 



27 



2(3 2 1 + (3 2 



(4) 



1 + /3 2 



(5) 



and (3 = g/m^q is the momentum transferred between the nucleons, is the 
charged-pion mass, and 7 will be discussed below. Details of this calculation (includ- 
ing graph-by-graph results), pion and nucleon cr-terms, and the renormalization-scale 
dependence of the coupling constants will be published elsewhere. 

This force is CD and remarkably simple in form, given the number of processes 
that contribute. Fourier transformation to configuration space yields 
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a 



4tt/2 Ltt 

I(z) = 2e-*[ln(s/2) + 7s - 7] + 2e z Ei(-2z) 



(ti • t 2 - t\ 4) a x -V Z (T 2 - V 2 



Hz) 



7 p * 

Ei(-z)(3 + -) + — 



(6) 
(7) 
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where Ei{—z) = — / 2 °° dt e~ l jt is the exponential integral, z = m n cr/h, 7^ is Euler's 
constant (0.577- • •), and the usual OPEP corresponding to charged-pion exchange is 
obtained by substituting e~ z for I(z) and dropping the factor of a/ it in brackets. Our 
renormalization prescription defines 7 and is discussed next. The first term in I(z) 
(in brackets) determines the asymptotic form, while the volume integral of I(z) is 
proportional to (| — 27). A positive value of 7 (~ 1) weakens both. 



The usual definition ||13|| of the ttNN coupling constant requires an extrapolation to 
the pion pole (i.e., q 2 = m 2 ) in the unphysical region of the NN scattering amplitude, 
with the residue defining the coupling constant and the difference between q 2 = 
and q 2 = m 2 defining the GT discrepancy [0]. This pole corresponds (since q° = 
in the static limit) to (3 2 = — 1. Although V nc (P) has a simple pole at (3 2 = — 1, the 
first term in V ni {(5) does not (the residue diverges) because of the logarithm induced 
by the infrared structure of the photon loops. This is reflected in the configuration 
space term: e~ z ln(z). Any constant multiple of e~ z in I(z) (i.e., a multiple of OPEP) 
can be arbitrarily transferred to OPEP with an appropriately redefined f 2 , since the 
sum of V W1 and V nc remains unchanged. Alternatively, the logarithmic (asymptotic) 
term in V n ^(f3) can be made to vanish at any convenient point. Straightforward 
development of V ny ((3) leads to 7 = (an even simpler result!). We choose, however, 
to remove the 7^ term in V nl (r) by performing a further finite renormalization and 
fixing 7 = 7#; this defines f 2 in the presence of EM corrections, and is analogous to 
the MS (Modified Minimal Subtraction) renormalization commonly used in Standard 
Model calculations. The bracketed (logarithmic) term in I(z) now vanishes at r = 
2h/m 7T c = 2.8 fm. As we have previously discussed, this weakens both the tail of the 
7T-7 potential and the volume integral of I(z) (by a factor of three). 



It is difficult to compare our result with previous calculations! 19, 20, 21 . None 



of the CD calculations^, |24|, |25], |26], [27] were complete, and few were written 
in an easily interpretable form. Some numerical results are available, however^]. 
The gauge invariance of our final result gives us confidence in its form. Figure (lg) 
is easily calculated in old-fashioned perturbation theory in the static limit (see Ref. 
|p8|), and a comparison checks our overall sign and factors. As an additional check, 



the Breit interaction |29| can be calculated using our conventions and yields the usual 
result. 



The tail of V nl (for r > 1.4 fm) was incorporated into the Nijmegen PWA of 
only the np data|30| on an equal footing with OPEP. A total of 4107 np data were 
fit, and the results are shown in Table I. Each entry lists the x 2 °f the fit and the 
fitted value of / c 2 (xl000). The first entry is the 1997 Nijmegen np PWA result [0. 
Simply adding V ni increases \ 2 by a factor of 8. This can be greatly reduced by 
refitting only f 2 , which is shown in the interim (next) entry. Refitting both f 2 and 
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Table 1: Nijmegen PWA fits. The vr ± iViV coupling constant, / c 2 (xl000), and the 
corresponding \ 2 °f fit for a variety of cases are indicated in the top row. The 
potentials are those included in the tail of the np force, while the quantities in brackets 
were those fit to produce the results listed directly below. "All" denotes that the 
parameterized interior region of the force was fit in addition to / 2 . 



Fit 


K[all] K + K 7 [/ c 2 ] Kr + K 7 [all] 


10 3 / c 2 

x 2 


74.96(34) 75.22 74.98(33) 
4223.6 4236.5 4222.8 



the phenomenological interior region (r < 1.4 fm) leads to the rightmost entry. Only 
those entries labelled "all" should be compared. 

Although it is necessary to increase slightly the strength of the tail of OPEP to 
compensate for the addition of V nl in the overall fit, there is negligible change in 
/ 2 , and only a tiny decrease in \ 2 ■ The np l S$ scattering length and effective range 
also show negligible change. We conclude that the addition of V ni has not affected 
any of the former conclusions of the Nijmegen PWA with regards to / 2 , quality of 
fit, or low-energy scattering observables. We emphasize that this does not mean 
that V nl is everywhere negligible, since that part with r < 1.4 fm is subsumed in 
the phenomenological interior region of the Nijmegen PWA. If a force were to be 
constructed (for all r) with an explicit Try component, the effect of the latter would 
be considerably larger. This np force is repulsive in S-waves for r < 3.7 fm, but 
attractive otherwise. A rough estimate using ad hoc form factors and several pp 
potentials shows that the magnitude of the np 1 Sq scattering length decreases by ~ 
0.67 fm. This is twice the strength of the effect found earlier from the double-seagull 



mechanism 1 27]. The deuteron energy is raised by roughly 60 keV. 



As a check of our procedures, we have also fitted the cases 7 = 7b ±1. These alter- 
natives to our renormalization scheme either add or remove a specified fraction (2a / tt) 
of OPEP from V^ 7 , and this must be compensated in the final fit by a corresponding 
change in / 2 . This is indeed found. Choosing 7 = (no finite renormalization) lowers 
1000/ 2 by 0.18, which is slightly more than half the quoted uncertainty. We elect, 
however, to use our preferred convention (7 = 7b), since that leaves / 2 essentially 
unchanged. 

In summary, we have calculated 19 Feynman graphs that produce the leading- 
order one-loop EM corrections to OPEP. The results are gauge invariant and remark- 
ably simple, and generate a static charge-dependent 7r-7-exchange potential. This 
potential, which has a nominal strength (a/7r)V^, was one of the few isospin-violating 
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mechanisms not previously incorporated into the tail of the NN force by the Nijmegen 
PWA. When incorporated using our renormalization prescription, there is negligible 
change in the low-energy parameters and the n ± NN coupling constant, f%, and 
only a tiny improvement in the quality of fit. 
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